(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     11153,        310]
NotebookOptionsPosition[      9795,        280]
NotebookOutlinePosition[     10169,        296]
CellTagsIndexPosition[     10126,        293]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["Closed Form Social Welfare according to HSV (2017) [EQ 30]", "Text",
 CellChangeTimes->{{3.806319800945203*^9, 
  3.806319831513025*^9}},ExpressionUUID->"17c1792e-9fc6-4459-819e-\
94a46c3d2ef8"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"W", "[", "\[Tau]_", "]"}], ":=", 
   RowBox[{
    RowBox[{"Log", "[", 
     RowBox[{"1", "-", "g"}], "]"}], "+", 
    RowBox[{"\[Chi]", "*", 
     RowBox[{"Log", "[", "g", "]"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "\[Chi]"}], ")"}], "*", 
     FractionBox[
      RowBox[{"Log", "[", 
       RowBox[{"1", "-", "\[Tau]"}], "]"}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", 
         FractionBox[
          RowBox[{"\[Sigma]", "+", "\[Tau]"}], 
          RowBox[{"1", "-", "\[Tau]"}]]}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"1", "-", "\[Tau]"}], ")"}]}]]}], "-", 
    FractionBox["1", 
     RowBox[{"1", "+", 
      FractionBox[
       RowBox[{"\[Sigma]", "+", "\[Tau]"}], 
       RowBox[{"1", "-", "\[Tau]"}]]}]], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "\[Chi]"}], ")"}], "*", 
     FractionBox[
      RowBox[{"Log", "[", 
       RowBox[{"1", "-", "\[Tau]"}], "]"}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "\[Psi]"}], ")"}], "*", 
       RowBox[{"(", 
        RowBox[{"\[Theta]", "-", "1"}], ")"}]}]], "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Psi]", "*", 
        RowBox[{"Log", "[", 
         RowBox[{"1", "-", "\[Tau]"}], "]"}]}], "+", 
       RowBox[{"Log", "[", 
        RowBox[{
         FractionBox["1", 
          RowBox[{"\[Eta]", "*", 
           SuperscriptBox["\[Theta]", "\[Psi]"]}]], "*", 
         SuperscriptBox[
          RowBox[{"(", 
           FractionBox["\[Theta]", 
            RowBox[{"\[Theta]", "-", "1"}]], ")"}], 
          RowBox[{"\[Theta]", "*", 
           RowBox[{"(", 
            RowBox[{"1", "+", "\[Psi]"}], ")"}]}]]}], "]"}]}], ")"}]}], "-", 
    RowBox[{
     FractionBox["\[Psi]", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", "\[Psi]"}], ")"}], "*", "\[Theta]"}]], "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "\[Tau]"}], ")"}], "-", 
       RowBox[{
        FractionBox[
         RowBox[{"\[Beta]", "*", "\[Delta]"}], "\[Gamma]"], "*", 
        FractionBox[
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Gamma]"}], ")"}], 
         RowBox[{"(", 
          RowBox[{"1", "-", 
           RowBox[{"\[Beta]", "*", "\[Delta]"}]}], ")"}]], 
        RowBox[{"(", 
         RowBox[{"1", "-", "\[CapitalTau]"}], ")"}]}]}], ")"}]}], "-", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Log", "[", 
        RowBox[{"1", "-", 
         FractionBox[
          RowBox[{"1", "-", "\[Tau]"}], "\[Theta]"]}], "]"}]}], "-", 
      FractionBox[
       RowBox[{"1", "-", "\[Tau]"}], "\[Theta]"]}], ")"}], "-", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "-", "\[Tau]"}], ")"}], "2"], "*", 
     FractionBox[
      RowBox[{"\[Upsilon]", "[", "\[Phi]", "]"}], "2"]}], "-", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "\[Tau]"}], ")"}], "*", 
       RowBox[{"(", 
        FractionBox[
         RowBox[{"\[Beta]", "*", "\[Delta]"}], 
         RowBox[{"(", 
          RowBox[{"\[Gamma]", "-", 
           RowBox[{"\[Beta]", "*", "\[Delta]"}]}], ")"}]], ")"}], "*", 
       FractionBox[
        RowBox[{"\[Upsilon]", "[", "\[Omega]", "]"}], "2"]}], "-", 
      RowBox[{"Log", "[", 
       FractionBox[
        RowBox[{"1", "-", 
         RowBox[{"\[Delta]", "*", 
          RowBox[{"Exp", "[", 
           RowBox[{
            FractionBox[
             RowBox[{
              RowBox[{"-", "\[Tau]"}], "*", 
              RowBox[{"(", 
               RowBox[{"1", "-", "\[Tau]"}], ")"}]}], "2"], "*", 
            RowBox[{"\[Upsilon]", "[", "\[Omega]", "]"}]}], "]"}]}]}], 
        RowBox[{"1", "-", "\[Delta]"}]], "]"}]}], ")"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "+", "\[Chi]"}], ")"}], "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        FractionBox["1", 
         FractionBox[
          RowBox[{"\[Sigma]", "+", "\[Tau]"}], 
          RowBox[{"1", "-", "\[Tau]"}]]], 
        RowBox[{"\[Upsilon]", "[", "\[Epsilon]", "]"}]}], "-", 
       RowBox[{"\[Sigma]", "*", 
        FractionBox["1", 
         SuperscriptBox[
          RowBox[{"(", 
           FractionBox[
            RowBox[{"\[Sigma]", "+", "\[Tau]"}], 
            RowBox[{"1", "-", "\[Tau]"}]], ")"}], "2"]], "*", 
        FractionBox[
         RowBox[{"\[Upsilon]", "[", "\[Epsilon]", "]"}], "2"]}]}], 
      ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.80631983550476*^9, 3.806319873985075*^9}, {
  3.806320162386161*^9, 3.806320530100754*^9}, {3.8063218709984827`*^9, 
  3.8063218830711346`*^9}, {3.8063223409874425`*^9, 3.8063228970208116`*^9}, {
  3.8063231052653255`*^9, 3.806323125048057*^9}, {3.806382861453251*^9, 
  3.8063828778340235`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"98e1c431-915d-4dce-882d-5635868d900a"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"g", "=", 
   FractionBox["\[Chi]", 
    RowBox[{"1", "+", "\[Chi]"}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalTau]", "=", "0.174"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Delta]", "=", "0.971"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Chi]", "=", "0.233"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Sigma]", "=", "2.00"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Upsilon]", "[", "\[Phi]", "]"}], "=", "0.036"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Upsilon]", "[", "\[Epsilon]", "]"}], "=", "0.164"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Upsilon]", "[", "\[Omega]", "]"}], "=", "0.003"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Theta]", "=", "3.124"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{" ", 
  RowBox[{
   RowBox[{
    OverscriptBox["\[Sigma]", "^"], "=", 
    FractionBox[
     RowBox[{"\[Sigma]", "+", "\[Tau]"}], 
     RowBox[{"1", "-", "\[Tau]"}]]}], ";", 
   RowBox[{"\[Psi]", "=", "0.65"}], ";", 
   RowBox[{"\[Gamma]", "=", "0.95"}], ";", 
   RowBox[{"\[Beta]", "=", "0.95"}], ";", 
   RowBox[{"\[Eta]", "=", "1"}], ";"}]}]}], "Input",
 CellChangeTimes->{{3.8063231158789606`*^9, 3.806323153652998*^9}, {
   3.8063232250923815`*^9, 3.806323241509183*^9}, {3.806323272948226*^9, 
   3.8063233861003895`*^9}, {3.8063805746726313`*^9, 3.806380617771017*^9}, {
   3.8063807909448285`*^9, 3.8063808169031415`*^9}, {3.8063812822652645`*^9, 
   3.8063812885122867`*^9}, {3.8063819980008073`*^9, 
   3.8063820239285326`*^9}, {3.806382477928646*^9, 3.806382552440867*^9}, 
   3.8063825983692093`*^9, {3.8063829310910273`*^9, 3.8063829449376087`*^9}, {
   3.8063830879299746`*^9, 3.806383093369239*^9}, {3.8063837588847275`*^9, 
   3.8063837617478137`*^9}, {3.8063841472048016`*^9, 3.806384147946565*^9}, {
   3.806384184348948*^9, 3.806384211347332*^9}, {3.8063843101228094`*^9, 
   3.8063843118432603`*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"f893670b-ec14-4f14-ba4e-38ae10d89da8"],

Cell["Calc the values", "Text",
 CellChangeTimes->{{3.8063835653537145`*^9, 
  3.8063835689695673`*^9}},ExpressionUUID->"af2077f4-a49f-495e-888d-\
ee4f151ae06a"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"W", "[", "0.181", "]"}]], "Input",
 CellChangeTimes->{{3.8063826296892147`*^9, 3.806382631288704*^9}, {
  3.8063835432447166`*^9, 3.8063835558107367`*^9}},
 CellLabel->"In[62]:=",ExpressionUUID->"33203c79-fe5f-44d1-a6d2-6879ccff7135"],

Cell[BoxData[
 RowBox[{"-", "1.1096142534016131`"}]], "Output",
 CellChangeTimes->{
  3.806382632879196*^9, {3.8063832619744315`*^9, 3.8063832764468126`*^9}, 
   3.8063837671330953`*^9, 3.806384318586108*^9},
 CellLabel->"Out[62]=",ExpressionUUID->"badbd759-0457-4b03-8a89-58febd229003"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"W", "[", "0.084", "]"}]], "Input",
 CellChangeTimes->{{3.8063835292181606`*^9, 3.80638353480189*^9}},
 CellLabel->"In[63]:=",ExpressionUUID->"87a04a85-b10e-45de-9f12-14c86b896201"],

Cell[BoxData[
 RowBox[{"-", "1.0873843603231008`"}]], "Output",
 CellChangeTimes->{3.8063835794737425`*^9, 3.80638377006516*^9, 
  3.806384320529488*^9},
 CellLabel->"Out[63]=",ExpressionUUID->"375750af-626f-42f7-80ca-cbfdbbbc89d4"]
}, Open  ]],

Cell["Relative change", "Text",
 CellChangeTimes->{{3.8063835739214973`*^9, 
  3.8063835757699165`*^9}},ExpressionUUID->"2255c27c-85d9-43bd-88a9-\
938d6c3b6149"],

Cell[CellGroupData[{

Cell[BoxData[
 FractionBox[
  RowBox[{"(", 
   RowBox[{
    RowBox[{"W", "[", "0.084", "]"}], "-", 
    RowBox[{"W", "[", "0.181", "]"}]}], ")"}], 
  RowBox[{"W", "[", "0.181", "]"}]]], "Input",
 CellChangeTimes->{{3.806383701979888*^9, 3.8063837332589054`*^9}},
 CellLabel->"In[64]:=",ExpressionUUID->"439608f7-4cd1-4f70-9ad1-6fb123b2aa8e"],

Cell[BoxData[
 RowBox[{"-", "0.020033892868954014`"}]], "Output",
 CellChangeTimes->{{3.806383723640975*^9, 3.806383734124629*^9}, 
   3.806383772718313*^9, 3.806384324638177*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"59f32a3a-39cd-46a1-bcb3-081a96c7b1bc"]
}, Open  ]]
},
WindowSize->{2048, 1109},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>1.4 Inherited,
FrontEndVersion->"11.3 for Microsoft Windows (64-bit) (March 6, 2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 200, 3, 48, "Text",ExpressionUUID->"17c1792e-9fc6-4459-819e-94a46c3d2ef8"],
Cell[761, 25, 4894, 143, 176, "Input",ExpressionUUID->"98e1c431-915d-4dce-882d-5635868d900a"],
Cell[5658, 170, 2105, 49, 341, "Input",ExpressionUUID->"f893670b-ec14-4f14-ba4e-38ae10d89da8"],
Cell[7766, 221, 161, 3, 48, "Text",ExpressionUUID->"af2077f4-a49f-495e-888d-ee4f151ae06a"],
Cell[CellGroupData[{
Cell[7952, 228, 258, 4, 40, "Input",ExpressionUUID->"33203c79-fe5f-44d1-a6d2-6879ccff7135"],
Cell[8213, 234, 287, 5, 46, "Output",ExpressionUUID->"badbd759-0457-4b03-8a89-58febd229003"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8537, 244, 204, 3, 40, "Input",ExpressionUUID->"87a04a85-b10e-45de-9f12-14c86b896201"],
Cell[8744, 249, 232, 4, 46, "Output",ExpressionUUID->"375750af-626f-42f7-80ca-cbfdbbbc89d4"]
}, Open  ]],
Cell[8991, 256, 161, 3, 48, "Text",ExpressionUUID->"2255c27c-85d9-43bd-88a9-938d6c3b6149"],
Cell[CellGroupData[{
Cell[9177, 263, 341, 8, 69, "Input",ExpressionUUID->"439608f7-4cd1-4f70-9ad1-6fb123b2aa8e"],
Cell[9521, 273, 258, 4, 46, "Output",ExpressionUUID->"59f32a3a-39cd-46a1-bcb3-081a96c7b1bc"]
}, Open  ]]
}
]
*)

